m 



In preparation for the Astro-physical Journal 

Three dimensional geometries and the analysis of Hn regions 

Kenneth Wood 1,2 , J. E. Barnes 1 , Barbara Ercolano 3 , L. M. Haffner 2 , R. J. Reynolds 2 , J. Dale 4 



^^ ■ kw25@st-andrews.ac.uk 

u 

On 1 abstract 

< 

in 

< 

o 



We compare emission line intensities from photoionization models of smooth and 
fractal shell geometries for low density Hn regions, with particular focus on the low- 
ionization diagnostic diagram [Nil] /Ha vs Ha. Building on previously published models 
and observations of Barnard's Loop, we show that the observed range of intensities and 
variations in the line intensity ratios may be reproduced with a three dimensional shell 
geometry. Our models adopt solar abundances throughout the model nebula, in contrast 
with previous one dimensional modeling which suggested the variations in line intensity 
ratios could only be reproduced if the heavy element abundances were increased by a 
factor of ~ 1.4. For spatially resolved Hn regions, the multiple sightlines that pierce and 
sample different ionization and temperature conditions within smooth and fractal shells 
produce a range of line intensities that are easily overlooked if only the total integrated 
intensities from the entire nebula model are computed. Our conclusion is that inference 
of Hn region properties, such as elemental abundances, via photoionization models of 
one dimensional geometries must be treated with caution and further tested through 
^^ , three dimensional modeling. 
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1. Introduction 

Analysis of emission lines from Hn regions provides information on the ionization state and tem- 
perature of the gas, elemental abundances within the nebula, and the Lyman continuum spectrum 
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of the ionizing source(s). Observations of Hn regions are usually interpreted using one-dimensional 
models where an external source ionizes a slab or an internal source ionizes a sphere or spherical 
shell and the glob al averages of the emi ssion line intensities and their ratios are compared with 
observations (e.g.. iPequignot et al.ll200ll ). Although employing one-dimensional radiation trans- 
fer, such models contain a wealth of very detailed physics describing atomic, molecular, and dust 
processes. 

Recent code development s have extended photoionization modeling to include th ree dimen- 
sional radiation transfer (e.g., lOch et alJll99a : Ercolano et al.l 120031 ; IWood et al.l 12004 ) and thus 
the ab ility to study the p rojected intensities from highly asymmetric nebulae. For example, mod- 
els by IWood et al.1 (J2005h of the Ha (A = 6563 A) and [Nil] (A = 6583 A) emission lines from 
the C Oph Hn region employing a fractal density distribution demonstrated the variety of line 
ratios that could be obtained. Ultimately the data and models in t his analysis wer e comp ared 
by taking radia l aver ages of the Ha intensity and [Nil] /Ha line ratio. lErcolano et al.l (|20071 ) and 
Ercolano et al.l (120101 ) explored the effects of a 3D distribution of ionizing sources on abundance 
determinations. These authors found that a centrally concentrated distribution of ionizing stars 
results in a higher effective ionization parameter in the nebula, compared to a more spread out 
distribution of the same stellar population. This can give errors of up to an order of magnitude in 
abundance, when strong line methods are used. However they found that abundance determina- 
tions based on direct measurements of temperature (by use o f nebular to auroral li ne diagnostics) 
are not affected. Models of star forming regions presented by lErcolano et al.l (|201ll ) incorporating 
3D hydrodynamics and photoionization further demonstrate the variety of line ratios present in 
spatially resolved 2D images. 

There are many uncertainties in photoionization modeling including atomic data, ionizing 
spectra, elemental abundances, and, the focus of this paper, the geometry of the gas. The two 
most popular m odels to describe Hii region data are that of a Stromgren sphere of gas ionized by 
a central source (|Stromgrenlll939l ) and the filling f actor model where ionized b lobs are surrounded 
by vacuum (e.g., see discussion and examples in lOsterbrock fc Ferlandl 120061 ) . Neither of these 
approaches represents the reality of an Hn region containing ionized and neutral gas with a range 
of temperatures, densities, and ionization states. Smooth shell models do have a range of temper- 
atures and ionization states that change as a function of distance from the ionizing source, but 
do not have neutral inclusions and ionized-neutral interfaces that naturally arise due to shielding 
in 3D geometries. Moving to 3D geometries gives a greater variety of temperature and ioniza- 
tion states within a m odel Hn region and a correspondingly larger spread in projected intensities 
(|Ercolano et alJboilh . 



Recently , usin g the CLOUDY photoionization code (as described in iFerland et al.1 Il998l ). 
O'Dell et al.l (|201ll ) analyzed new and archival emission line data from Barnard's Loop and con- 
cluded that there may be an increase in elemental abundances. Their models assumed that the 
abundances of all elements were increased in unison and as such did not take into account differ- 
ent processes producing the elements. In this paper we construct spherical shell and 3D fractal 



shell geometries for Barnard's Loop that demonstrate the observations may be reproduced without 
appealing to abundance variations. 



2. 3D photoionization simulations 

The ph o toion ization models are constructed using the 3D Monte Carlo code described by 



Wood et al.l fj2004j ). The code discretizes the Hn region density on a 3D linear Cartesian grid 
and computes the electron temperature and ionization structure of H, He, C, N, O, Ne, and S, 
the elements that dominate the cooling in low density Hn regions. We do not consider photons 
with energies above 54 eV, the ionization potential of He + , which is appropriate for the ionizing 
sources and consistent with observations which demonstrate there is ve ry little ionized helium in 
low density Hn regions and diffuse ionized gas (JRevnolds &: Tuftdll995l ). Due to grid resolution, 
we do not compute emissivities for cells where the neutral fraction is greater than 0.25. Such cells 
represent ionized-neutral interfaces and our current 3D Cartesian grid cannot spatially resolve the 
rapid change in ionization state at the interface. 



We follow lOPell et al.l (|201ll ) and start off by constructing spherically symmetric models 
comprising a central point source illuminating a shell of uniform density, n = 3cm~ 3 , with inner 
and outer radii R{ n = 46 pc and Rr , u t = 5 5 pc. We fix the ionizing luminosity at Q = 2.5 x 



10 49 s _1 as estimated bv lOPell et al.l ( 201ll ) for the stars believed to be ionizing Barnard's Loop. 
Photoionization models are constructed for source temperatures 29000 K and 32000 K, spanning 
the estimated temperature range of the io nizing stars. The ioni zing spectra for our simulations 
were taken from the WM-BASIC library of ISternberg et al.l (J2003J) assuming solar abundances and 
log g = 4. For 3D models we use the same inner and outer r a dii an d convert the uniform density 
shell to a fractal g eometry using the algorithm o f ElmegreenJ (119971 ) as described in several of our 
recent papers (e.g. JMathis et al.ll2002l : IWood et al.ll2005l ). This algorithm leaves a fraction, / sm0 oth, 
in a smooth, uniform density component, and redistributes the remainder into hierarchical clumps. 
We adopt / smoo th = 1/3, so the lowest density in our fractal Barnard's Loop models is n = 1cm -3 . 



Images of Barnard's Loop (e.g., Figure 1 in lODell et al.ll201ll ) suggest a shell thickness of 
around lOpc. This is indeed achieved in our simulations because the smooth models are ionization 
bounded at around 54pc. For the fractal models, their smooth component is ionized out to the 
outer edge at 55pc, but because of the dominance of ionized higher density clumps the overall 
appearance is that of a thin fragmented shell (see Figure 1 below). 

In all the simulations presented below the elemental abundances are He/H= 0.1, C/H= 1.4 x 
10" 4 , N/H= 6.5 x 10~ 5 , 0/H= 4.3 x 10~ 4 , Ne/H= 1.17 x 10~ 4 , and S/H= 1.4 x 10~ 5 . Our 
photoionization models provide emissivities and hence projected intensities for various emission 
lines and for this paper we compute maps of Ha (A = 6563 A) and [Nil] (A = 6583 A). Our code 
also provides the intensi ty of the two closely spaced lines of [Sn] (A = 6716 A and A = 6731 A). 
However, as discussed by lO Dell et al.1 (1201 ll ). we can shift the model [Sii]/Ha intensity ratio by an 



unknown amount because we do not know the dielectronic recombination rates for sulfur. Therefore, 
in this paper we focus our attention on the [Nil] /Ha vs Ha diagram. 



3. Results 

Figure 1 shows the Ha intensity and [Nil] /Ha line ratio maps for the 29000 K source for the 
smooth and fractal shell models (the corresponding maps for the 32000 K source are qualitatively 
similar). We can immediately see the Ha limb brightening for the smooth shell and the more 
fragmented, but still ring-like structure for the fractal shell model. Figures 2 and 3 present the 
intensity and line ratio maps as scatter plots for the diagnostic diagram [Nil] /Ha vs Ha. Clearly 
a wide range of values occur for both Ha intensity and line ratios as we now explain. 

The ionizing luminosity, shell density and radial extent, are such that the smooth shell is 
ionization bounded (i.e., all ionizing photons are absorbed before the outer edge of the shell), 
resulting in very high temperatures at the edge of the ionized volume. This is shown in Figure 4 
and arises because only the highest energy photons can penetrate to the outermost edges of the 
ionized volume, hence producing the highest temperatures in the nebula. In contrast the fractal 
shell has a smooth density component which is n = 1 cm" 3 and this component is ionized, unless 
shadowed by dense clumps. Therefore the temperatures in the smooth component of the fractal 
model are smaller than for the corresponding radius in the smooth shell model (see Figure 4). 
Hence there is not the very large [N n]/Ha associated with low Ha intensities as seen at the edges 
of the smooth shell. 

The two ionizing source temperatures for the smooth model yield line ratios that nicely bracket 
the observations. The double valued nature of the [Nil] /Ha vs Ha diagram for the smooth model 
shown in Figure 2 is explained with reference to Figure 4 which shows intensity cuts and also the 
radial temperature profile within the shell. Towards the edge of the shell the Ha intensity rises and 
falls with radius due to limb brightening. However, the [Nil] /Ha ratio increases with radius because 
sightlines towards the edge of the shell are sampling higher and higher temperatures. Hence the 
double- valued nature of this diagnostic diagram. The points with low Ha and large [Nil] /Ha arise 
at large radius while the low Ha and [Nil] /Ha are towards the inner regions of the shell. It is 
unlikely that such double-valued diagnostic diagrams will be observed because of the difficulty of 
accurately subtracting foreground and background emission from the warm ionized medium, the 
low density ionize d component of the interstellar medium (see discussion of background effects in 



Wood et all 12003 ) . 



The right hand panels in Figures 2 and 3 show all impact parameters between 46 pc < R < 
54 pc, a range appropriate for comparison to Barnards loop. It is clear from fi gure 2 that a smooth 



density structure is unable to match the observations, a result also found by lO'Dell et al.1 (1201 ll ) 
which led them to suggest abundance variations as a method for better modelling the region. 
However figure 3 shows that if a fractal 3D model is used, it is possible to match the observations 



without the need to alter elemental abundances. 

The fractal shell models presented in Figures 1 and 3 display a wider range of Ha intensities 
due to the density contrasts within the simulation, with densities in the fractal shell often exceeding 
n > 30cm~ 3 . In addition, the values of [Nil] /Ha extend to lower values than the uniform density 
models. This corresponds to the cells in the simulation with the lowest densities (n = 1 cm -3 ) being 
ionized and the corresponding lower temperatures in the smooth ionized component as described 
above. Figure 4 shows the range of temperatures in the fractal models, nicely explaining the lower 
[Nil] /Ha compared to the smooth models. The [Nil] /Ha vs Ha diagram for the 29000 K source 
(red dots in Figure 3) shows that for all but the lowest H a intensities the line ratio lies in the range 



0.2 < [Nil] /Ha < 0.25, bracketing the data presented by lO 'Dell et all (|201ll ) 



Our 3D models do include the high temperatures that occur at ionized-neutral interfaces (see 
Figure 4). However our Cartesian grid does not adequately spatially resolve these boundaries in 
the fractal models. In our simulations the percentage of ionized grid cells we ignore with neutral 
fractions above 0.25 is less than 10% of the ionized volume. We anticipate that, while important, 
better resolution of interfaces within the simulation will not change our overall conclusions regarding 
the effects of 3D geometries on the line ratio diagnostic diagrams of Hn regions studied in this paper. 
The resolution effects of our Cartesian grid do mean that we cannot make accurate predictions of, 
for example, the [Oi] (A = 6300 A) emission line which is temperature sensitive and through charge 
exchange is tied to the ionization state of hydrogen. The [Oi]/Ha intensity ratio is larger in diffuse 
ionized gas than in Hn regions, suggesting that high temperat ures possibly associa ted with interface 



emission is important in widespread diffuse ionized gas (e.g.. lHausen et al.ll2002l ). 



The results we have presented for analytically produced fractal densities are also seen in pho- 
toi onization simulations of hydrodynamical simulations. The star forming region models presented 



by lErcolano et al.1 (120111 ) clearly demonstrate the same effects of a wide range of observed line ratios 
in spatially resolved images. See their figures 4 and 5 for a very nice example of diagnostic diagrams 
where the synthetic line ratios span a range which in cludes values that are usually associated with 



shock ionization. However, the lErcolano et al.1 (|201ll ) results are from purely photoionization simu- 



lations, again demonstrating that ID analysis can lead to misunderstandings and misclassifications 
of the physical processes producing observed emission line ratios. 



4. Summary 

Through 3D photoionization modeling we have demonstrated that the observed line ratios 
from Barnard's Loop may be reproduced without appealing to increased elemental abundances. 
We have focussed on the [Nil] /Ha vs Ha diagnostic diagram since accurate predictions of [Sn] 
emission is not possible due to unknown dielectronic recombination rates for sulfur. The simple 
explanation of our result is that 3D geometries provide a range of temperatures and ionization 
states and hence a range of emission line intensities different from the total integrated intensity 



from ID shell models (see also lErcolano et al.l l201ll ). However, the ID shell models presented 



in Figure 2 do in fact demonstrate that it is possible to reproduce the observations with smooth 
models when rays piercing different impact parameters are considered. A spatially re s olved analysis 



of the projected intensities from smooth models was not presented by lO'Dell et al.l (|201ll ) — this 
is the main point at which our analyses differ and the reason for our different conclusions regarding 
elemental abundances within Barnard's Loop. The uncertainties introduced by different atomic 
databases and 3D geometries make it difficult to determine absolute Hn region abundances from 
analysis of diagnostic diagrams. However, it will be possible to determine relative abundances and 
abundance gradients from photoionization modeling of Hll regions at different spatial locations 
within a galaxy. 

We thank Bob O'Dell and Gary Ferland for constructive criticism of our manuscript 
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Fig. 1. — Ha intensity maps (left) and [Nil] /Ha line ratio maps (right) for the smooth shell (upper 
panels) and fractal shell (lower panels) models. The axes are labeled in units of parsecs. The values 
of the intensities and line ratios are shown as scatter plots in Figures 2 and 3 
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Fig. 2. — Diagnostic diagrams for the smooth shell models. Red and blue dots are results for 
the 29000 K and 32000 K sources respectively. In the left panel, each dot represents a different 
sightline through the model nebula, corresponding to the individual pixels in the two dimensional 
maps presented in Figure 1. The right panel shows the models for impact parameters in the range 
46 pc < R < 54 pc, corresponding to the observations of Barnard's Loop (diamonds). The double- 
valued pattern in the left hand panel arises from the increasing temperature with radius in the 
shell and the corresponding rise then fall of the Ha intensity with radius. The highest values for 
the line ratios occur at the outer edge of the nebula (see Figure 4). 
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Fig. 3. — As for Figure 2, but for the fractal shell models. Left panel shows intensities and 
line ratios for all impact parameters R < 55 pc, right panel restricts impact parameters to be in 
the range 46 pc < R < 54 pc. Notice the wide range of Ha intensities that arise from the range of 
densities, temperatures, and ionization fractions sampled by different sightlines through the models. 
The [Nil] /Ha values do not reach the high values as for the smooth models because the smooth 
component of the fractal models is almost fully ionized and therefore at a lower temperature than 
the corresponding radius in the smooth models. 
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Fig. 4. — Left panels: intensity cut across the center of the Ha (solid line) and [Nil] /Ha (dashed line) 
line ratio maps from Figure 1 for the smooth (upper) and fractal (lower) models. The intensities 
and line ratios have been normalized to their maximum values. This explains the double-values 
obtained in the Ha vs [Nil] /Ha diagnostic diagram for the smooth models in Figure 2. Upper right 
panel: radial temperature profile of the smooth model showing the rise in temperature towards the 
edge of the Hn region. Lower right panel shows the temperature at all radial points in the fractal 
shell model. The concentration of points between 5500 K and 6000 K corresponds to the smooth 
component of the fractal shell — notice these temperatures are lower at large radii than for the 
smooth shell model. Although the shell physically extends to i? ut = 55pc, the ionizing luminosity 
is such that for the smooth shell the edge of the Hn region occurs around r ~ 54pc, hence the 
rapid rise in temperature as the gas turns neutral and only highest energy photons can reach the 
outermost locations. The fractal shell's smooth component is ionized out to the outer edge of the 
shell (low temperatures) and the high temperatures present in the lower right panel arise at the 
ionized-neutral interfaces that occur throughout the fractal shell. 



